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Resumen 


Este trabajo presenta el proceso necesario de adecuación de datos de precipitación 
obtenidos de imágenes de radar meteorológico mediante SIG para su empleo en 
el modelado hidrológico distribuido. El objetivo principal ha sido la evaluación 
de la calidad de la información de radar a través de las simulaciones del 
modelo hidrológico y la comparación con la simulación generada con datos de 
pluviómetros. Para llevar a cabo lo anterior, se simularon tres eventos. Se incluyó 
la avenida extrema que se presentó en la cuenca del río Escondido, México, en abril 
del 2004; se utilizaron datos horarios del radar meteorológico ubicado en la base 
aérea Laughling Texas, Estados Unidos, que consisten en archivos de tipo binario 
en un listado de valores de precipitación en mm. La calidad de los datos de radar 
se llevó a cabo mediante un análisis de correlación con datos de pluviógrafos de 
estaciones automatizadas. Los resultados obtenidos indican que en la comparación 
radar-pluviógrafo existen pares de datos en los que mientras en el radar se estima 
un valor “cero”, en el pluviógrafo se está registrando un valor mayor que cero 
o viceversa. Por otra parte, la varianza explicada de los datos de precipitación 
estimados a partir de información de radar, con respecto a las mediciones de los 
pluviógrafos en las estaciones consideradas en el análisis, fue menor a una escala 
horaria que en la escala diaria. Finalmente, las simulaciones realizadas con datos de 
precipitación-radar para eventos extremos muestran que el hidrograma reproduce 
acertadamente el hidrograma observado. En cambio, en la simulación efectuada 
con datos pluviométricos, tanto el pico como la forma del hidrograma, no fueron 
simulados de manera adecuada. 


Palabras clave: NEXRAD, radar meteorológico, modelo distribuido, SIG, 
CEQUEAL. 


Introducción 


Hoy en día, los hidrólogos y especialistas invo- 
lucrados en la solución de diversos problemas 
de recursos hídricos y ambientales reconocen la 
importancia de incorporar una alta resolución 
a la variabilidad espacial y temporal de la 


precipitación como uno de los factores más 
importantes en las aplicaciones hidrológicas 
(Pathirana et al., 2003). Adicionalmente, en el 
modelado hidrológico existe un creciente interés 
en el uso de datos de precipitación estimados 
a partir de imágenes de radar. En efecto, la 
estimación de la precipitación mediante datos 
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generados por radar es particularmente útil 
en zonas montañosas, donde la topografía 
accidentada y la altitud del terreno limitan 
significativamente la operación de las redes 
climatológicas debido a dificultades de acceso 
(Hossain et al., 2004). Sin embargo, también 
es precisamente en zonas montañosas donde 
se llegan a presentar dificultades para la 
estimación de la precipitación con imágenes 
de radares, ya que la señal emitida por el radar 
puede verse interferida por las montañas, 
limitando el análisis a la zona no interferida 
por obstáculos físicos. 

Actualmente, las estimaciones de preci- 
pitación a partir de información de radar 
se utilizan en el modelado del fenómeno 
lluvia-escorrentía, lo cual permite maximizar 
la aplicación de modelos hidrológicos y es- 
pecialmente los de tipo distribuido. Las 
aplicaciones incluyen el pronóstico hidrológi- 
co, la determinación de zonas de inundación y 
la evaluación de los recursos hídricos (Durrans 
et al., 2002; Jessen et al., 2005; Krajewski et 
al., 2006). Además de la aplicación directa de 
los radares en meteorología, la aplicación de 
éstos a fenómenos hidrológicos requiere de 
una calibración, la cual consiste en obtener 
relaciones de reflectividad para diferentes 
intensidades de precipitación, específicamente 
dentro de la zona de influencia (Rosengaus, 
1995). Este proceso se debe realizar a través de 
una red de pluviógrafos con una resolución 
temporal del orden de minutos. 

Con respecto a la modelación hidrológica, 
Neary et al. (2004) utilizaron datos de radar 
y datos de pluviógrafo con el modelo 
hidrológico HEC-HMS para determinar la 
aplicabilidad de la información de radar con 
modelos no distribuidos o globales. Estos 
autores encontraron que los escurrimientos 
simulados con datos procedentes de radar, en 
general fueron subestimados por el modelo 
hidrológico y los gastos picos observados 
no fueron simulados adecuadamente. Con 
respecto a la aplicación de datos de radar en 
modelos hidrológicos distribuidos, Kalin y 
Hantush (2006) utilizaron datos de precipita- 


ción estimados a través de radar, así como datos 
de pluviógrafos, conjuntamente con el modelo 
SWAT (Neitsch et al., 2002). Particularmente se 
interesaron en la exploración del potencial de 
los datos de radar como una fuente alterna de 
los datos de precipitación obtenidos a través de 
pluviógrafos. 

Por otro lado, cabe recordar que el 
principio de funcionamiento del radar meteo- 
rológico consiste en la emisión de pulsos 
electromagnéticos hacia la atmósfera, y los 
equipos detectan los ecos o reflejos de energía 
de los conglomerados de lluvia presentes en la 
atmósfera. De esta manera se pueden ubicar 
tormentas caracterizando su estructura, así 
como cuantificar la intensidad de la lluvia. La 
cobertura de un radar meteorológico típico 
es de aproximadamente 240 km de radio. Sin 
embargo, la cobertura real es menor debido 
a la curvatura de la Tierra y a la orografía 
del lugar donde se encuentre instalado. Los 
radares no miden directamente el volumen 
de agua, sino la reflectividad Z (ecuación (1)) 
o energía reflejada por las gotas de agua. Las 
principales fuentes de incertidumbres en los 
datos arrojados por el radar meteorológico 
son debidas precisamente a que el radar no 
mide directamente la precipitación, la cual 
es convertida en intensidad R (ecuación 
(2)) mediante la relación Z-R (ecuación (3)), 
deducida a partir de un análisis estadístico 
entre los datos del radar y datos de estaciones 
convencionales (p.e. pluviómetros). Otra fuen- 
te de incertidumbre se debe a que la precipita- 
ción a nivel de la superficie es deducida a 
través de medidas del radar a ciertas alturas 
por encima del suelo, y la variabilidad de la 
precipitación en la altitud puede causar una 
sobrestimación o subestimación (Jayawick- 
reme y Hyndman, 2007). Esto implica que para 
utilizar los datos de precipitación inferidos 
a partir de imágenes de radar es necesario 
estimar la incertidumbre antes de utilizar la 
información en el modelado hidrológico. 
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Donde: 


N(D) es una función de distribución del 
número de gotas de agua. 

V(D) es una función de distribución del 

volumen de las gotas de agua. 

D es el diámetro de las gotas de agua. 

es el máximo diámetro detectado de 


gotas de agua. 


máx 


En la ecuación (3), los parámetros a y b 
dependen de la calibración del radar y se han 
obtenido mediante un análisis de correlación 
utilizando datos de radar, así como datos 
de precipitación obtenidos a través de 
pluviógrafos. Por ejemplo, en el cuadro 1 se 
presentan algunos valores de los parámetros a 
y b utilizados para transformar la reflectividad 
en precipitación (Bedient et al., 2000; Javier et 
al., 2007). 

El presente estudio muestra el proceso 
necesario de adecuación de datos de preci- 
pitación obtenidos de imágenes radar 
meteorológico mediante un SIG (sistema de 
información geográfica) para su empleo en el 
modelado hidrológico distribuido. El objetivo 
principal ha sido la evaluación de la informa- 
ción de radar a través de las simulaciones del 
modelo hidrológico y la comparación con 
la simulación obtenida con base en datos de 


pluviometría. Particularmente, se simularon 


tres eventos en la cuenca del río Escondido, 
México: dos ordinarios en 2003 y la avenida 
extrema que se presentó en abril del 2004. 


Metodología 
Zona de estudio 


La cuenca del río Escondido se localiza al norte 
de México y se origina en la Sierra del Burro; por 
la margen izquierda recibe al río San Antonio, 
para posteriormente cruzar por la ciudad de 
Piedras Negras antes de unirse por la margen 
derecha con el río Bravo. La cuenca tiene un 
área de 3 175 km? hasta la estación hidrométrica 
Villa de Fuentes. El tipo de vegetación en la 
cuenca es básicamente matorral submontano. 
La corriente principal nace en una altitud de 
1 000 metros sobre el nivel del mar (msnm) 
y recorre una longitud de 160 km hasta una 
elevación de 245 msnm, donde se localiza la 
estación hidrométrica. La precipitación media 
anual en la cuenca varía de 400 a 600 mm. 

La zona de estudio cuenta con tres estacio- 
nes meteorológicas convencionales y un obser- 
vatorio meteorológico (figura 1), los cuales son 
operados por la Comisión Nacional del Agua 
(Conagua) y la información está disponible en 
la base de datos CLICOM (Climatic Computing 
Project) a través del Servicio Meteorológico 
Nacional (SMN). Adicionalmente existe una 
red de estaciones automatizadas sobre la 
margen derecha del río Bravo que son operadas 
por la Comisión Internacional de Límites y 
Aguas entre México y Estados Unidos (CILA). 
Con respecto a la información hidrométrica, el 


Cuadro 1. Relaciones Z-R utilizadas para la estimación de la precipitación con base en imágenes 


de radar meteorológico. 


z Parámetro 
Ubicación del radar , Relación Z-R 
a 
Llanura 300 1.4 R =300Z'* 
Montaña 119 (0177 e = MEZR 
Tropical 250 1.2 R=250Z!'? 


Fuentes: Bedient et al., 2000; Javier et al., 2007. 
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registro de caudales medios diarios medidos 
en la estación hidrométrica Villa de Fuentes 
data de 1932 a la fecha y se obtuvo de http: / / 
www.ibwc.state.gov. Los caudales se obtienen 
a partir de las curvas de aforo de la estación, 
la cual cuenta con un vertedor triangular para 
la medición del flujo en época de estiaje y una 
canastilla para aforar con molinete hidráulico 
en crecidas. El promedio de los caudales 
diarios máximos de 1932 a 2007 es de 59.3 m*s”. 
Sin embargo, durante el periodo de registro se 
han presentado varios eventos extraordinarios 
(cuadro 2), siendo el máximo registrado el 
ocurrido en abril de 2004. 


Datos NEXRAD 


En México existe una red de radares meteo- 
rológicos equipados con tecnología Doppler y 
operados por la Conagua a través del SMN. Las 
estaciones cuentan con sistema ininterrumpible 
de energía, sistema de protección contra 


Cuadro 2. Caudales máximos registrados en el río 
Escondido de 1932 a 2007. 


Año Caudal máximo (m? s”) 
1936 190 
1957 218 
1964 371 
1969 193 
2004 883 


Fuente: http: / /www.ibwc.state.gov. 


incendios y sistema de comunicaciones vía 
satélite con el centro colector de datos que se 
ubica en las instalaciones del SMN en la ciudad 
de México. Actualmente, la cobertura del 
territorio nacional no se encuentra a un 100%, 
por lo que la cuenca del río Escondido no está 
cubierta por el radio de influencia del radar 
mexicano más cercano. 

Debido a lo anterior se utilizaron datos 
de un radar meteorológico localizado en la 
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Figura 1. Zona de estudio. 
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base aérea Laughling Texas, Estados Unidos. 
Este radar es del tipo WSR-88D (Weather 
Surveillance Radar 1988 Doppler), conocido como 
la generación futura de radares meteorológicos 
o NEXRAD (next-generation weather radar). La 
información generada por este tipo de radares 
fue diseñada específicamente para cubrir 
la parte continental de Estados Unidos con 
base en una proyección polar estereográfica 
(Reed y Maidment, 1999), conocida como 
HRAP (Hydrologic Rainfall Analysis Project). 
Con respecto a la disponibilidad de los datos 
fuente utilizados, éstos únicamente abarcan 
parte de la frontera norte de México y alcanzan 
a cubrir la totalidad de la cuenca del río 
Escondido (figura 2). La información se obtuvo 
de la NOAA (National Oceanic and Atmospheric 
Administration) de los Estados Unidos. Se 
utilizaron archivos horarios del 1? de enero de 
2002 al 31 de diciembre de 2004 (37 224 archivos 
y 239 faltantes). Esta información se encuentra 
en archivos horarios de tipo binario, en un 
listado de valores de precipitación en mm. El 
formato de los datos es de tipo real de cuatro 
octetos (con siete cifras significativas) para 
cada celda de la rejilla HRAP. Para la zona de 


Figura 2. Imagen de precipitación obtenida con el radar 


sobre la cuenca del río Escondido a las 8:00 am del 4 de 
abril del 2004. 


estudio se extrajo la información contenida en 
coordenadas HRAP de las columnas 489 a 517 
y renglones 110 a 127. 

Para determinar la calidad de la informa- 
ción se llevó a cabo un análisis estadístico, 
el cual consistió en comparar mediante un 
análisis de regresión lineal los datos de radar 
contra los datos de precipitación de tres 
estaciones automatizadas  (pluviógrafos), 
las cuales se localizan sobre la margen 
derecha del río Bravo y son operadas por la 
CILA. En el cuadro 3 se muestran tanto las 
estaciones convencionales como las estaciones 
automatizadas utilizadas. En la combinación 
de datos de precipitación se utilizaron pares de 
datos con la precipitación estimada a partir de 
las imágenes de radar de la malla HRAP que 
coincidía con la ubicación de cada uno de los 
pluviógrafos y para cada uno de los intervalos 
analizados. En el análisis se emplearon datos 
de precipitación a intervalos de tiempo 
horarios y diarios, verificando que los datos 
de radar y de pluviógrafo estuvieran en escala 
de tiempo UTC. Cabe mencionar que sólo 
se consideraron pares de datos que tuvieran 
registros de precipitación concomitantes. 

Finalmente, para hacer coincidir la malla 
de puntos de precipitación de radar con la 
malla de cuadros del modelo hidrológico 
utilizado, se realizó un proceso de remuestreo 
del sistema HRAP utilizando una función 
de interpolación lineal. Esto consistió en 
estimar la esperanza matemática de los 
valores horarios de los puntos precipitación 
de radar circunscritos en cada cuadro en 
que se discretizó la cuenca. Se automatizó el 
proceso de remuestreo, así como la obtención 
de archivos de precipitación en el formato que 
requería el modelo hidrológico. Lo anterior se 
efectuó mediante el diseño e implementación 
de un módulo hidrogeomático inserto en el 
SIG Idrisi (Quentin et al., 2007). 


Modelo hidrológico 


La evaluación de la calidad de los datos de 
radar se basó tanto en el análisis estadístico 
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Cuadro 3. Estaciones terrenas convencionales y automáticas utilizadas en el análisis. 


Estación Ubicación 
Registro de datos 
Nombre Tipo Longitud Latitud 
i 1 
Zaragoza da -100.929 28.429 Diario 
(pluviómetro) 
Convencional Ea 
Allende E ES -100.859 28.354 Diario 
(pluviómetro) 
Convencional E 
Allende S a -100.833 28.333 Diario 
(pluviómetro) 
El Moral aa -100.634 28.898 Horario 
(pluviógrafo) 
Aut tizad. 
Piedras Negras aia -100.506 28.713 Horario 
(pluviógrafo) 
Río Escondido Ada -100.557 28.668 Horario 
(pluviógrafo) 


correspondiente, así como en la precisión de 
los resultados obtenidos en las simulaciones 
realizadas con el modelo hidrológico. 
El modelo 
distribuidos utilizado en esta investigación 
es el modelo CEQUEAU (Morin ef al., 1998). 


Este modelo es de tipo determinístico y de 


hidrológico de parámetros 


simulación continua, que discretiza la cuenca 
en estudio en una malla de cuadros. Considera 
las variaciones espaciotemporales tanto de las 
características físicas como de la precipitación 
(Bá et al., 2001). La simulación hidrológica 
del modelo se basa en tres recipientes: el 
recipiente lagos y ciénagas, el recipiente 
suelo y el recipiente subsuelo. Consta de dos 
partes que definen la forma en que se obtienen 
el flujo o caudal a la salida de la cuenca. La 
primera parte es la función de producción que 
modela el movimiento vertical del agua en los 
recipientes. Para cada cuadro que discretiza 
la cuenca, el modelo asigna la precipitación 
obtenida del 
respecto a las distancia, de las tres estaciones 


promedio, ponderado con 
más cercanas a dicho cuadro. La segunda 
parte, definida como función de transferencia, 
es la encargada de modelar el movimiento del 
agua de elemento en elemento o de cuadro por 
cuadro. El coeficiente de transferencia de cada 


cuadro está relacionado con las características 
hidráulicas parciales predominantes del flujo, 
es decir, la capacidad de amortiguar la onda 
mediante el depósito en la red. Un buen 
índice de la capacidad de amortiguación en 
una parte del cuadro es el área de agua libre 
en el cuadro. Un gran lago, por ejemplo, 
absorbe el flujo, pero su influencia también 
depende del área de la cuenca aguas arriba. El 
modelo CEQUEAU utiliza los coeficientes de 
transferencia para cada cuadro, considerando 
el área cubierta por el agua y la zona aguas 
arriba de dicho cuadro. 

La validación de las simulaciones se 
llevó a cabo con base en el criterio numérico 
propuesto por Nash y Sutcliffe (1970). Este 
criterio cuantifica la precisión de los resultados, 
relacionando los caudales observados Qo y 
los caudales calculados Qc expresados por 
la ecuación (4). El dominio matemático de 
existencia del valor del coeficiente de Nash es 
de -oo hasta 1, donde la unidad representa la 
simulación perfecta: 


Y (Qc;- o)? 

NTD=1-% —-o<NTD<1 (4) 
Y (o; - Lo)” 
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La implementación del modelo hidroló- 
gico requirió de: a) un archivo de datos fisio- 
gráficos, b) un archivo de datos hidro- 
meteorológicos y c) un archivo con los 
parámetros del modelo. Para crear el archivo 
fisiográfico se utilizó un Modelo Digital de 
Elevación (MDE) generado por interferometría 
de radar, como un nuevo módulo 
hidrogeomático que trabaja en el SIG Idrisi, 


así 


y desarrollado para automatizar el proceso 
de extracción de los sentidos del flujo y la 
ocupación del suelo de cada cuadro de la malla 
que discretiza la cuenca en formato CEQUEAU 
(Guerra-Cobián, 2007). Se construyeron los 
archivos de datos hidrometeorológicos: el 
primero con los datos de los pluviómetros 
localizados en la zona de estudio, el segundo 
con los datos del radar meteorológico en el 
formato del modelo hidrológico. Por último, 
se efectuó la calibración y validación con base 
en el ajuste de los parámetros del modelo 
hidrológico considerando dos escenarios. 
En el primer escenario sólo se utilizaron 
los datos de las estaciones convencionales 
(datos de pluviómetro). El segundo escenario 
de simulación se realizó con los datos 
de precipitación estimada con radar. El 
procedimiento de calibración se realizó en dos 
etapas: en la calibración primaria se utilizó 
la técnica de prueba y error, y en la segunda 
etapa, o calibración óptima, se empleó la 
herramienta de optimización de parámetros 
contenida en el modelo CEQUEAU. Se 
calibraron los parámetros del modelo que 
controlan el flujo vertical del agua (alturas de 
vaciado, umbrales de agua en los recipientes, 
etcétera), así como los coeficientes del flujo de 
los recipientes. 


Resultados y discusión 

Análisis estadístico de los datos de radar 

El proceso de verificación de la calidad de 
los datos de radar se llevó a cabo mediante 


datos 
correspondientes de pluviógrafos, los cuales 


un análisis de contraste con los 


se muestran en la figura 3. Con la dispersión 
existente de los datos se aprecia que la 
precipitación estimada por el radar difiere 
significativamente de la precipitación medida 
por los pluviógrafos para los intervalos de 
tiempo analizados. Estos resultados coin- 
cidieron en los tres pluviógrafos seleccionados 
y de acuerdo con Neary et al. (2004), esta 
dispersión no es atribuible totalmente al 
radar, sino que probablemente se debe a las 
fuentes de incertidumbre en las mediciones 
de pluviógrafos. En efecto, esta incertidumbre 
puede atribuirse a un error sistemático 
en la medición de viento y evaporación, 
así como un inadecuado funcionamiento 
mecánico. Otra fuente de incertidumbre es 
probablemente ocasionada por la diferencia 
entre las áreas de toma de datos, ya que la 
precipitación obtenida por medio del radar 
está promediada en una celda cuadrada de 
16 km?” y el dato del pluviógrafo es puntual. 
En el cuadro 4 se presentan los resultados del 
análisis de correlación entre los datos de radar 
y de pluviógrafo. Se observa que los valores 
del coeficiente de determinación (R?) de las 
tres estaciones consideradas en el análisis, a 
intervalo de tiempo horario, son menores que 
con intervalo de tiempo diario (cada 24 horas). 
De las estaciones analizadas, la estación río 
Escondido fue la que presentó el mejor ajuste 
bajo unintervalohorario(R?=0.60), y la estación 
Piedras Negras con intervalo de tiempo diario 
(e = 
intervalo diario coinciden con los resultados 
encontrados por Kalin y Hantush (2006) en 
análisis similares, ya que reportan valores de 


0.92). Los resultados obtenidos con 


R? entre 0.87 y 0.91 para el mismo intervalo de 
tiempo. Asimismo, los resultados del presente 
estudio indican que, en la comparación 
radar-pluviógrafo, existen pares de datos 
concomitantes que, mientras en el sensor del 
radar se estima un valor “cero”, en el sensor 
del pluviógrafo se está registrando un valor 
mayor que cero o viceversa. Esto podría estar 
relacionado directamente con la capacidad del 
sensor en detectar la precipitación. Además, 
contrario a los resultados obtenidos en otras 
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Figura 3. Dispersión de datos de precipitación horaria y diaria radar-pluviógrafo. 
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Cuadro 4. Resultados del análisis estadístico. 


Estación Intervalo de tiempo R? Modelo de ajuste 
Horario 0.32 Y = 0.36x + 1.35 
El Moral 
Diario 0.78 Y = Més 175 
Horario 0.60 Y =0.722x + 1.13 
Río Escondido 
Diario 0.91 Y =0.93x + 1.44 
Horario 0.44 Y = 0.73x + 0.92 
Piedras Negras 
Diario 0.92 Y =1.11x - 0.149 
Cuadro 5. Resultados de la simulación hidrológica. 
Caudal medio Caudal máximo Coeficiente de Nash 
(mi s”) (més”) (adimensional) 
Observado 0.98 20.80 - 
Iculad 
can 8.91 5.96 -0.272 
Evento 1 pluviómetros 
Calado isa 2.937 13.37 0.356 
precipitación-radar 
Observado 0.58 9.06 - 
Eaeaco on 13.61 6.63 -0.500 
Evento 2 pluviómetros 
o 1.04 6.70 0.132 
precipitación-radar 
Observado 10.03 883.0 - 
coco 5.19 61.6 0.079 
Evento 3 pluviómetros 
Ellis 8.35 884.5 0.974 
precipitación-radar 


cuencas por Jayawickreme y Hyndman (2007), 
la estimación de la precipitación con radar en 
los píxeles que coinciden con los pluviógrafos 
fue sobrestimada. 


Proceso de remuestreo de los datos de radar 


La función de remuestreo del software Idrisi 
fue empleada para obtener una restauración 
geométrica de las imágenes. Se utilizó para 
georregistrar las imágenes en un sistema 
de referencia mediante las coordenadas de 
una serie de puntos de control en un archivo 
existente, así como en el nuevo sistema de 
referencia deseado, convirtiendo el archivo 


creado al nuevo sistema de referencia por 
medio de una función cartográfica polinomial 
lineal. Los resultados gráficos del proceso 
utilizado se muestran en la figura 4 y 
consisten en una imagen de puntos del sistema 
HRAP remuestreada en formato raster. Este 
procedimiento permitió analizar gráficamente 
la distribución espacial de la precipitación 
en la cuenca del río Escondido. Igualmente, 
se observa que entre las 7:00 y 8:00 horas 
del día 4 de abril del 2004, la distribución 
de la precipitación no cubre la totalidad de 
la superficie de la cuenca. Esta precipitación 
no fue registrada en las estaciones de terreno 
debido a que éstas se encontraban alejadas 
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Figura 4. Remuestreo de datos de precipitación HRAP a precipitación CEQUEAU. 


del centro y área de influencia de la tormenta. 
De acuerdo con los datos de las estaciones 
terrenas, se registró un valor máximo de 15 
mm de precipitación, mientras que el valor 
máximo estimado por el radar fue de 101.1 
mm (tomando como referencia la precipitación 
acumulada entre las 7:00 y 8:00 am del 4 de 
abril del 2004). 


Simulación hidrológica 


La simulación se efectuó en el periodo com- 
prendido de 2002 a 2004, lapso en el que se 
cuenta con la información del radar. Además, 
la calibración se llevó a cabo utilizando los 
datos de precipitación estimada a partir de 
imágenes de radar. Lo anterior se debió a que 
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sólo se cuenta con tres pluviómetros, y éstos no 
están distribuidos de forma adecuada dentro 
de la cuenca. Ello ocasionó que no se pudiera 
efectuar la calibración del modelo únicamente 
con datos de pluviómetros. 

Con respecto a los resultados obtenidos 
de la simulación, la figura 5 muestra la 
comparación de los hidrogramas observado y 
simulado en la cuenca del río Escondido, así 
como sus respectivos hietogramas. Para los 
eventos 1 y 2 se observa que los hidrogramas 
no fueron simulados satisfactoriamente. Sin 
embargo, en el evento 1, la simulación con 
datos de precipitación-radar es conside- 
rablemente mejor que la simulación con datos 
de pluviómetros. En cambio, en el evento 
3, el hidrograma obtenido en la simulación 
realizada con datos de precipitación-radar 
reproduce adecuadamente el hidrograma 
observado. También se observa en el hidro- 
grama obtenido en la simulación efectuada 
con datos de pluviómetros que no fue posible 
reconstruir de manera adecuada el caudal que 
se presentó durante el evento. Los valores del 
coeficiente de Nash, así como los caudales 
medios del año 2004 y máximo anual se 
muestran en el cuadro 5. El coeficiente de Nash 
que se obtuvo en el proceso de calibración 
corresponde a la avenida máxima del mes de 
abril del 2004 (evento 3). Con respecto a los 
escurrimientos durante los años 2002 y 2003, 
fueron prácticamente cero, excepto los eventos 
1 y 2. Lo anterior es una característica de la 
cuenca del río Escondido, ya que se encuentra 
en una zona semidesértica, donde las lluvias 
son de tipo convectivo, asociadas con frentes 
fríos y ocasionalmente con eventos ciclónicos. 
Además, el caudal base es prácticamente nulo 
para varios meses del año. 

Resulta evidente que la deficiencia del 
modelado con base en información pluvio- 
métrica se debió a que la distribución de 
la tormenta en la cuenca se localizó en una 
zona donde no se contaba con dispositivos de 
medición de precipitación. 

Por otro lado, cabe mencionar que a 
pesar de que el modelo CEQUEAU no está 


diseñado para trabajar con un formato de 
datos de precipitación estimada a partir de 
radar meteorológico, gracias a que el modelo 
discretiza la cuenca en una malla de cuadros, 
ésta se adecuó fácilmente a la malla de los 
datos de precipitación estimados por imágenes 
de radar. 


Conclusiones 


Las estimaciones de la precipitación con base 
en imágenes de radar son particularmente 
útiles en zonas donde las tormentas son aisla- 
das o de tipo convectivas y, específicamente, si 
la red de estaciones meteorológicas no cubre 
eficientemente estas zonas. Complementa- 
riamente, el uso de los sistemas de infor- 
mación geográfica (SIG) como parte del pre- 
procesamiento de la información requerida 
en el modelado hidrológico es hoy en día una 
herramienta indispensable. En efecto, los SIG 
son necesarios para el tratamiento y manejo 
de las imágenes de precipitación-radar debido 
a que permiten la visualización y edición de 
la información, así como la extracción de los 
datos de precipitación mediante el desarrollo 
de algoritmos computacionales insertos en el 
SIG. 

Los datos de precipitación-radar utilizados 
en el modelado hidrológico permitieron 
simular tres eventos. Entre éstos en el 
evento 3 (abril de 2004), considerado como 
extraordinario, fue simulado de forma 
ampliamente satisfactoria tanto el caudal pico 
de la avenida como la forma del hidrograma. 
En cambio, para los eventos ordinarios 1 y 
2 no fue posible reconstruirlos de manera 
adecuada. Sin embargo, en forma general, 
se considera una modelación apropiada 
para eventos potencialmente generadores 
de inundación con daños asociados. Los 
resultados obtenidos mostraron que a nivel 
horario existe una baja correlación entre 
los datos de precipitación medida por 
pluviómetro y la estimada por imágenes de 
radar. En cambio, a nivel diario se obtuvieron 


buenas correlaciones. Finalmente, es de suma 
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Figura 5. Resultados de la simulación hidrológica efectuada en CEQUEAU con datos de pluviómetro y RADAR. 


importancia en el modelado hidrológico 
conocer las capacidades y limitaciones de los 
datos de precipitación estimados a través de 
radar, ya que esto permitirá su aplicación de 
una manera más eficiente y facilitará la toma 
de decisiones oportunas. 
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Abstract 


GUERRA-COBIÁN, V.H., BÁ, K.M., QUENTIN-JORET, E., DÍAZ-DELGADO, C. $ 
CÁRSTEANU, 4.A. Use of NEXRAD information for hydrological modeling in watersheds 
with sparse rain gauge networks. Water Technology and Sciences, formerly Hydraulic 
engineering in Mexico (in Spanish). Vol. IL, No. 3, July-September, 2011, pp. 35-48. 


This paper presents the process to adjust precipitation data obtained from weather radar images 
using GIS, for use in distributed hydrologic modeling. The main objective was to evaluate the 
quality of radar data using hydrologic model simulations and comparisons with simulations 
generated from rain gauge data. To this end, three events were simulated, including the 
flood that occurred in the Escondido River basin, Mexico, in April 2004. Hourly data were 
used from weather radar located in Laughling AFB, Texas, USA, which consist of binary 
files from a list of values of precipitation in mm. The quality of radar data was performed 
by a correlation analysis with data from automated rain gauge stations. The results indicate 
that data pairs exist in the radar-rain gage comparison, in which the radar is calculated as 
“zero” and the gauge registers a value greater than zero, or vice versa. Moreover, the explained 
variance between the precipitation data calculated from radar information and rain gauge 
measurements from the stations in the analysis was less at the hourly scale than the daily 
scale. Finally, simulations with radar precipitation data for extreme events show that the 
simulated hydrograph accurately reproduces the observed hydrograph. On the other hand, 
in the simulated hydrograph for rainfall data, neither the peak discharge nor the hydrograph 
shape was adequately simulated. 


Keywords: NEXRAD, meteorological radar, distributed model, GIS, CEQUEAU. 
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